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Abstract. We apply the non-linear diffusive shock acceleration theory in order to describe the 
properties of two supernova remnants, SN 1572 (Tycho) and SN 1604 (Kepler). By analyzing 
the multi-wavelength spectra, we infer that both Tycho's and Kepler's forward shocks (FS) are 
accelerating protons up to ~ 500 TeV, channeling into cosmic rays more than 10 per cent of their 
kinetic energy. We find that the streaming instability induced by cosmic rays is consistent with the 
X-ray morphology of the remnants, indicating a very efficient magnetic field amplification (up to 
~ 300/i G). In the case of Tycho we explain the y-ray spectrum from the GeV up to the TeV band 
as due to pion decay produced in nuclear collisions by accelerated nuclei scattering against the 
background gas. On the other hand, due to the larger distance, the y-ray emission from Kepler is not 
detected, being below the sensitivity of the present detectors, but it should be detectable by CTA. 
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INTRODUCTION 

Modelling particle acceleration at SNRs is one of the main goal to explain the origin 
of Galactic cosmic rays (CRs) in the context of the so-called supernova paradigm [see 
e.g. 1]. This paradigm requires that SNRs are able to accelerate nuclei up to energies 
as high as a few times 10 6 GeV, converting a fraction ~ 10% of the SNR kinetic 
energy into CRs. There is an increasing amount of evidence that shocks in young SNRs 
can indeed reach the required efficiency. Unfortunately, all these evidence are indirect 
and when considered individually can be also explained by other mechanisms that do 
not require efficient acceleration. A clear example of such an ambiguous situation is 
represented by y-ray emission. In the last few years several SNRs have been detected 
both in the GeV and in the TeV band but, in spite of this increasing amount of data, 
the question whether this emission is due to hadronic (through the decay of neutral 
pions produced in nuclear interactions between accelerated nuclei and the background 
plasma) or leptonic processes (due to inverse Compton and/or bremsstrahlung) is still 
debated [see e.g. 2, for a general discussion on this topic]. A possible way to discriminate 
between the leptonic (i.e. inefficient) and the hadronic (i.e. efficient) scenario is studing 
the effects that particle acceleration produces at all observable wavelengths, rather that 
focus the attention only on the y-ray emission, simultaneously using multiple set of data 
to constrain the model. 

Here we apply the non-linear diffusive shock acceleration (NLDSA) theory in order 
to explain the non-thermal emission and morphological properties of two similar SNRs, 
namely Tycho and Kepler. Tycho, in particular, has been recently detected in y-rays by 



Fermi-LAT [3] and VERITAS [4], and can be considered one of the most promising 
object where to test the shock acceleration theory and hence the CR-SNR connection. 
We presented a detailed model apply ed to Tycho in [5]. Kepler, on the other hand, has 
not been detected in y-ray band, yet, but it is very similar to Tycho in many respects, and 
we show that the predicted y-ray spectrum should be detectable by CTA. 

COUPLING PARTICLE ACCELERATION WITH REMNANT 

EVOLUTION 

Both Tycho and Kepler have been established to be remnants of a type la SN, Tycho 
thanks to observation of the scattered-light echo [6] while Kepler based on the O/Fe 
ratio observed in the X-ray spectrum [7]. They have similar ages and similar values 
for the radio spectral index, i.e. 0.65 and 0.64, respectively. Interestingly also the ratio 
of power emitted at 10 keV and at 1 GHz, i.e. vF v (@10keV)/vF v (@lGHz), is very 
similar, being 56.9 for Tycho and 52.9 for Kepler. Finally, both remnants present very 
thin filaments in non-thermal X-ray emission. 

In order to make the model as simple as possible, in both cases we assume that the 
remnants expand into an uniform circumstellar medium (CSM) with proton number 
density no (which we leave as a free parameter) and temperature Tq = 10 4 K. We model 
the remnant evolution by following the analytic prescriptions given by [8] for Type la 
SNe, assuming a SN explosion energy E$n — 10 51 erg and one solar mass in the ejecta, 
whose structure function is taken as oc (v/v e j)~ 7 . The radial structure of density and 
temperature profiles is then calculated by assuming that the shocked CSM is roughly in 
pressure equilibrium. 

On top of this SNR evolution, the spectrum of accelerated particles is calculated ac- 
cording to the semi-analytic kinetic formalism put forward in [9] and references therein, 
which solves self-consistently the equations for conservation of mass, momentum and 
energy along with the diffusion-convection equation describing the transport of non- 
thermal particles for quasi-parallel, non-relativistic shocks. The injection is regulated by 
a free parameter, <^ n j , in such a way that all particles from the thermal plasma with mo- 
mentum p > %mjPth,2> with p t ^2 the typical thermal momentum downstream, start the 
acceleration process. 

A crucial role in our model is played by the magnetic field amplification induced by 
the super- Alfvenic streaming of relativistic particles upstream of the shock. We model 
this magnetic amplification as in [10]. The large magnetic fields predicted by the res- 
onant amplification have two main consequences: 1) the magnetic pressure upstream 
becomes comparable to, or even larger than, the thermal plasma pressure and, reduc- 
ing the compressibility of the plasma, affects the shock compression factor. 2) When 
the magnetic field is amplified the velocity of the scattering centers, which is generally 
neglected with respect to the shock speed, may be significantly enhanced [10]. When 
this occurs, the total compression factor felt by accelerated particles may be appreciably 
reduced and, in turn, the spectra of accelerated particles may be considerably softer. We 
explicitly include these effects assuming that the scattering centers moves with a speed 
equal to the Alfven speed in the amplified magnetic field. Once the magnetic field struc- 
ture is known, we can compute the spectrum of accelerated electrons at the shock, which 



is assumed to be proportional to the proton spectrum f e (p) = K ep f p (p), up to a maxi- 
mum momentum /?^ m ax determined by the synchrotron losses in the amplified magnetic 
field, where the spectrum presents a squared exponential cut-off oc exp[— (p/p^max) 2 ]- 
The evolution of the electron and proton spectrum downstream of the shock is computed 
taking into account adiabatic losses for protons and adiabatic and synchrotron losses for 
electrons. 

RESULTS. We consider the following radiative processes: 1) synchrotron emission 
of relativistic electrons; 2) thermal and non-thermal electron bremsstrahlung; 3) inverse 
Compton scattering (ICS) of electrons on CMB radiation, local IR dust emission and 
Galactic IR+optical light; 4) emission due to the decay of 71° produced in hadronic 
collisions. Because the explosion energy and the mass of the ejecta are fixed a priori, 
in order to fit the observed spectra and the remnant angular size we can only vary three 
parameters: the number density of the upstream medium, no, the injection efficiency of 
protons, £i n j, and the electron to proton normalization, K ep . 

In Fig. 1 we show our best fit of the photon spectrum produced by the superposition 
of all the radiative processes outlined above, comparing it with the existing data. The 
overall agreement is quite good. Quite interestingly, for both Tycho and Kepler our best- 
fitting returns no = 0.3 cm -3 and <^ n j =3.7 which implies that the total amount of kinetic 
energy channeled into CRs is ~ 12%. The inferred values for the circumstellar medium 
implies that Tycho and Kepler have a distance of 3.3 and 7.4 kpc from the Sun. 

Synchrotron radiation fits well the observed X-ray emission assuming an electron to 
proton ratio K ep = 1.6 x 10 -3 for Tycho and 2.8 x 10 -3 for Kepler. The projected X-ray 
emission profile is shown in Fig. 2, where it is compared with available Chandra data. 
In the case of Tycho the predicted profile shows a very good agreement with the data 
while for Kepler the model overpredicts the emission in the inner part of the rim. This 
could be due to a deviation of the shock geometry from the perfectly spherical symmetry 
assumed in our model. The sharp decrease of the emission behind the FS is due to the 
rapid synchrotron losses of the electrons in a magnetic field as large as ~ SOOjUG. 

In the y-ray band the decay of 71° produced in hadronic collision is the dominant pro- 
cess. In particular, we predict a slope for accelerated protons oc E~ 2 ' 2 . In the case of 
Tycho this slope well accounts for Fermi-LAT and VERITAS data within the experi- 
mental errors. The predicted proton spectrum shows a cut-off around /? max = 470 TeV/c. 

ICS of relativistic electrons cannot explain the observed y-ray emission for two 
different reasons. First, the strong magnetic field produced by the CR-induced streaming 
instability forces the number density of relativistic electrons to be too small to explain 
the y-ray emission as due to ICS on the ambient photons. Second, even if we arbitrarily 
reduced the magnetic field strength, enhancing at the same time the electron number 
density in order to fit the TeV y-rays with ICS emission, we could not account for 
the GeV y-rays because both the spectral slope and the flux would be incompatible 
with the Fermi-LAT data. Also non-thermal electron bremsstrahlung has to be ruled out, 
because it provides a flux two order of magnitudes lower than the Fermi-LAT detection, 
and cannot be arbitrarily enhanced without over-predicting both the TeV and the X- 
ray emission. Kepler, due to its larger distance, it is fainter than Tycho and therefore 
not detected in the y-ray band by the current generation of telescopes. However, the 
predicted TeV flux should be observable by the Cerenkov Telescope Array. 
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FIGURE 1. Spatially integrated spectral energy distribution of Tycho (left panel) and Kepler (right 
panel). The curves show synchrotron emission (thin dashed), thermal electron bremsstrahlung (dot- 
dashed), 71° decay (thick dashed) and ICS (dotted) as calculated within our model. The total emission 
is showed by the solid curve. The experimental data for Tycho are, respectively: radio from [11]; X- 
rays from Suzaku (courtesy of Tom Tamagawa), GeV y-rays from Fermi-LAT [3] and TeV y-rays from 
VERITAS [4]. For Kepler: radio from [11] and X-rays from [12]. 
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FIGURE 2. Projected X-ray emission at 1 keV for Tycho (left panel) and at 4-6 keV for Kepler (right 
panel) as a function of the distance from the shock position. The Chandra data points are from [13] 
(for Tycho) and from [14] (for Kepler). The dashed lines shows the projected radial profile of synchrotron 
emission while the solid lines are the same but convolved with the Chandra point spread function (~ 0.5"). 
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